Dynamic criticality in driven disordered systems: 
Role of depinning and driving rate in Barkhausen noise* 
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Jozef Stefan Institute, P.O. Box 3000, 1001 Ljubljana, Slovenia 

We study Barkhausen noise in a diluted two-dimensional Ising model with the extended domain wall 
and weak random fields occurring due to coarse graining. We report two types of scaling behavior 
corresponding to (a) low disorder regime where a single domain wall slips through a series of positions 
when the external field is increased, and (b) large disorder regime, which is characterized with 
nucleation of many domains. The effects of finite concentration of nonmagnetic ions and variable 
driving rate on the scaling exponents is discussed in both regimes. The universal scaling behavior 
at low disorder is shown to belong to a class of critical dynamic systems, which are described by a 
fixed point of the stochastic transport equation with self-consistent disorder correlations. 



I. INTRODUCTION 



Barkhausen signal is a physical phenomenon of large 
practical importance for noninvasive material character- 
ization technique and for magnetic recording. It con- 
sists of series of irregular pulses (noise) which are re- 
lated to reversal of the magnetization by slowly increas- 
ing (e.g., from negative saturation) external field in ferro- 
magnetic alloys and amorphous ferromagnets. Recently 
Barkhausen noise (BN) has been studied as an exam- 
ple of complex behavior in driven disordered systems far 
from equilibrium. It has been found that BN obeys 
a scale-free behavior without tuning of a parameter, a 
feature which is peculiar to the dynamic critical systems 
|^. In particular, most often measured are distributions 
of the area s of Barkhausen pulses, which obeys a power- 
law up to a cut-off according to D(s) ~ s~ Ts , the dura- 
tions of pulses t with distribution P(t) ~ t~ Tt , and the 
energy released in a pulse n with D(n) ~ nT Tn . The 
dissipated power spectrum was found to decay with fre- 
quency / as S(f) ~ f~^. A short summary of the ex- 
perimental results taken from the literature is given in 
Table 1. It should be noted that the scaling exponents 
obey (within error bars) the following scaling relations 
(t s - l)D s = (t„ - l)D n = n - 1, where D s and D n 
are the fractal dimensions associated with the size s and 
energy n of pulses, respectively. The measured scaling 
exponents can be grouped in (at least) three distinct uni- 
versality classes, which we may characterize according to 
the exponent t s as 1.33, 1.40, and 1.70, respectively. The 
effects of varying driving rate have been studied system- 
atically in certain compounds in Refs. [|]j7|,||,[l(| . It has 
been recognized that increased driving rate leads to in- 
creased occurrence of large pulses and thus to a contin- 



uous decrease of the scaling exponents. In addition, a 
careful analysis of the experimental results suggests that 
scaling exponents may depend on sample composition 
and annealing (compare, for instance, MC3 and Dl in 
Table 1), and on the exposing of a magnetostrictive sam- 
ple to tensile stress, which leads to different structure of 
domains 0, as opposed to unstressed sample ]I(J (see 
(DS) and (DB) in Table 1). 

Proximity of the mean-field universality class (corre- 
sponding to t s — 1.66, which is close to measured w 1.7 
in Table 1) was understood as being a consequence of 
dipolar forces and demagnetizing fields in d — 3 spatial 
dimensions jl0|. Recently occurrence of additional two 
universality classes (characterized with r s = 1.54 and 
t s = 1.30) |ll]] has been attributed to two-dimensional 
character of domain-wall motion, which was observed in 
a class of anisotropic soft magnets. Such systems are 
for instance ribbons made of Fe-B-Si jlJJ^] and Fe-Co- 
B compounds j7|,p|JTc|] , in which domains appear to have 
stripe structure with extended domain walls parallel to 
the magnetization. The two universality classes have 
been related |ll| to a single domain wall motion and to 
a multidomain structure, respectively, in a 2-dimensional 
prototype model with local random fields. Models with 
other types of disorder in d = 2 dimensions 13 -lq] have 
considered BN in multidomain structure. 

In the present paper we extend the study of Ref. M] 
to consider Barkhausen noise in two regimes correspond- 
ing to a single domain and multi-domain structure, re- 
spectively, in a diluted Ising model in d = 2 with vary- 
ing concentration of nonmagnetic ions c and fixed (weak) 
random fields. We use infinitely slow driving in order to 
study effects of the varying concentration c on the statis- 
tics of BN avalanches. We also discuss the effects of finite 



* Dedicated to the memory of my father 
Electronic address: Bosiljka.Tadic@ijs.si 



1 



driving rates in both regimes (see also [|14j). Finally, we 
show that the critical behavior in the regime with a sin- 
gle domain wall can be related to a class of driven critical 
systems which are described by the stochastic transport 
equation with correlated quenched currents fHI . 

In Section II we introduce the model and determine 
two disorder regimes. In Section III numerical results of 
distributions and their scaling exponents are presented 
for varying concentration c and driving rates. In Section 
IV we discuss the analogy between low-disorder regime 
and stochastic transport equation with quenched random 
currents, and in Section V we give a short summary and 
discussion of the results. 

TABLE I. Values of the scaling exponents measured in var- 
ious substances: (MCI) amorphous iron, (MC2) Ni95Al3Mn2, 
(MC3) FersBiaSig, all from Ref.l; (BG) VITRO VAC 6025 
X, Ref.6; (E4) Fe 3 oNi45B 15 Siio, Ref.4; (DS) Fe 64 Co 21 B 15 un- 
der tensile stress, Ref.8; (TR)1.8% SiFe for driving rate 0.15, 
Ref.5. (Dl) Fe 73 Coi 2 Bi5 for driving rate 0.1, Ref. 7; (DB) 
Fe 6 4Co 2 iBi 5 for driving rate 0.16, Ref. 10. Except in (MCI) 
and (MC2) which used wires, in the rest of experiments thin 
ribbon or strip samples were used. 
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II. MODEL AND BARKHAUSEN AVALANCHES 

We consider an Ising model on the square lattice in two- 
dimensions assuming that local random fields hi are gen- 
erated by coarse-graining from an original disorder in the 
presence of the external magnetic field H |lj] : 

« = - E J U S * S i - »< + " :S '< ■ (D 

<i,j> i 

Jij — 1 is a constant interaction between nearst-neighbor 
spins Si — ±1, and J^.j = if at least one of the spins is 
zero. A randomly distributed fraction c of sites represents 
nonmagnetic ions with Si — 0. A Gaussian distribution 
of hi is assumed with zero mean and width /. We assume 
a small random field variance / < 0.3 (see Ref. @ for 
effects of varying random fields), and vary concentration 
of nonmagnetic ions c. 

Following Ref. flllfl , we consider two disorder regions: For 
low values of concentrations c < c*(/) we consider ex- 
tended domain wall of flipped spins along < 11 > di- 
rection on the square lattice (linear size and the wall is 
L). Periodic boundaries are applied in the direction of 
the wall, and the boundary opposite to the wall is left 
free. Here c*(f) is determined as the largest disorder at 



which domain-wall depinning still occurs. For high dis- 
order c > c*(/) the domain wall remains pinned for all 
values of the driving field, and many domains of reversed 
spins are nucleated inside the system. In this region of 
disorder periodic boundary conditions in both directions 
can be used equally well as the boundary conditions de- 
scribed above, leading to the same avalanche statistics. 
In Fig. 1 we show two examples of emerging cluster struc- 
ture of flipped spins in low-disorder and high-disorder 
regime. The cluster structure at low disorder reflects the 
anisotropy due to initial conditions. The anisotropy is 
still important at high disorder, however, the occurrence 
of many domains that block each other spatial extent 
becomes dominant feature that determines the scaling 
properties of Barkhausen noise in this region. 
We first apply infinitely slow driving (i.e., increments of 
the driving field are adjusted to the minimum local field 
in the system). It should be stressed that these driving 
conditions arc useful for the purely theoretical purposes, 
as opposed to finite driving rates that are applied in real 
experiments. An example of the simulated BN signal for 
a half of the hysteresis loop is given in Fig. 2. The corre- 
sponding time series of the values of the field adjustments 
is given in the inset to Fig. 2. It is clear from Fig. 2 that 
after a certain transient time the consecutive values of 
the weakest local field in the system vary very little and, 
at the same time, the main contribution to the avalanche 
statistics comes from exactly that part of the hysteresis 
loop. 

III. CONCENTRATION AND DRIVING-RATE 
DEPENDENCE OF SCALING EXPONENTS 

In this Section we study distributions of avalanche sizes 
and durations by keeping fixed weak random field / <0.3 
at each occupied site and vary concentration of nonmag- 
netic ions c. We first apply infinitely slow driving in or- 
der to differentiate effects of varying disorder from those 
of varying driving rate. The concentration dependence 
for finite driving rate in high-disorder region was studied 
in Ref. [jlij. At low concentrations c the concept of a 
single domain wall persists, although a small fraction of 
spins may flip in the interior of the system near vacancies. 
The domain wall moves in a series of slips and eventu- 
ally depinning occurs at critical driving field H c (c). In 
the applied driving conditions the depinning may occur 
for concentrations c < c*(/), whereas for larger disor- 
der the built-in domain wall remains pinned for all val- 
ues of the driving field. For c > c*(/) the domains of 
flipped spins are nucleated inside the system and their 
walls move at different instances of time. In this way 
subsequent avalanches block each others extent, leading 
to long-range correlations in space and to memory effects. 
The two regions with different behavior are schematically 
shown in the phase diagram in Fig. 3. 
In the high disorder region HDR above c*(f) we measure 
the distribution of avalanche durations, which is shown 
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in Fig. 4 for varying c and fixed / = 0.001. Both the 
slope and cutoff of the distributions depend on c. In 
the inset are shown the extracted scaling exponents for 
slopes of the distribution of durations (r t ), and sizes (t s ), 
which arc found to increase with c. By diluting the spin 
lattice clusters of connected spins become smaller and 
more ramified, so that large avalanches are less probable 
by fixed other conditions (i.e., fixed / and driving rate). 
This leads to the increase of slopes of the avalanche dis- 
tributions, as shown in Fig. 4. 

Similarly, we find increase of the scaling exponents in the 
low disorder region (LDR) when c is varied in the inter- 
val < c < c*(/), where a single domain wall exists. 
Here we have measured the exponent of avalanche sizes 
t s (c), and derived Tt(c) from the scaling relations that 
are discussed in detail in Ref. |TI|. The exponents are 
also shown in the inset to Fig. 4. It is interesting to note 
that in this region of disorder the avalanches are self- 
affine (cf. Fig. 1). We measure the anisotropy exponent 
C, which is defined as < £j_ >^ ~ 1$, where i\\ and £± are 
avalanche extents in the direction parallel and perpendic- 
ular to the wall, respectively. For instance, for c = 0.05 
we find C = 1-23 ±0.09. 

It is expected that the depinning transition along the 
line H c (c) belongs to the same universality class as the 
depinning in c = limit for low random-field type of dis- 
order, which we studied in Ref. |D] ]. When the disorder 
is tuned to a critical value along the line c*(/) the phase 
transition occurs between low-disorder and high-disorder 
regime, which we expect to be in the same universality 
class as the one at critical random-field variance /* for 
c = (see Ref. II for details). 

By applying a finite driving rate, i.e., increasing the ex- 
ternal field always by a given constant amount AH in- 
dependently of values of local fields in the system, weak 
pinning centers become overdriven. Therefore avalanche 
sizes increase compared to the case of infinitely slow driv- 
ing both due to coalescence of different clusters and by 
faster propagation of a single cluster due to reduced pin- 
ning. As a result scaling exponents of avalanches decrease 
with increasing driving rate (see Ref. p5|). The effects of 
increased driving rate are similar in both LDR and HDR, 
however, the boundary between these regions in the 
phase diagram moves towards lower disorder (i.e., nucle- 
ation of new domains starts at lower disorder compared 
to the case of zero driving rate and interaction between 
domain walls prevents depinning transition). In the inset 
to Fig. 4 we have also shown values of the scaling expo- 
nents in HDR for driving rate 5h = AH/H max = 0.01, 
taken from Ref. [jTij. The concentration dependence of 
the exponents at finite driving rates becomes approxi- 
mately linear (cf. inset to Fig. 4 and Ref. |l4j). Another 
important feature of the scaling behavior by the finite 
driving rate is that the cutoffs of the distribution curves 
increase and obey the finite-size scaling [ pT[ . In this way a 
subcritical dynamic system for c*(/) < c < c p (f), where 
c p (f) is the percolation threshold line, becomes critical 



for a range of values of driving rates Jig ]. This is the 
reason for ubiquity of self-similar Barkhausen noise in 
highly disordered alloys, since the experimental setups 
always use finite frequency of the external field cycles. 
Dependence of the scaling exponents of the finite driving 
rates has been nicely demonstrated in Refs. [ ||7],|8|,[lO| in 
1.8% Si-Fe and Fe-Co-B compounds. It should be noted 
that both numerical jl4 15 and experimental observa- 
tions |||,^|J|Jl^] agree that for large driving rates 5h the 
exponents decay linearly with Sh. However, for small 
driving rates deviation from linear law occurs, making 
extrapolation to zero rate difficult. The dependence of 
the measured exponents in the literature on the concen- 
tration of nonmagnetic ions can be established qualita- 
tively. For instance, the duration exponent r t = 1.56 
measured in Fe73Coi2Bi5 (for driving rate 0.3 in reduced 
units) in Ref. ]7j is lower than r t — 1.82 measured in 
diluted Fe78B 13 Sig (for driving rate 1Hz) in Ref. Q, in 
a qualitative agreement with our predictions (cf.inset to 
Fig. 4). However, a more quantitative comparison would 
be only possible if all other conditions, i.e., driving rates 
and annealing of samples were kept under strict control. 



IV. UNIVERSALITY CLASS OF DRIVEN 
DISORDERED SYSTEMS 

In the low-disorder regime (/</*, c < c*(/)) apart 
from a weak c dependence the critical behavior at in- 
finitely slow driving can be identified by the exponents 
t s = 1.58, r t — 1.89 and £ = 1.23. These exponents ap- 
pear to be close to those obtained in the ricepile model 
with stochastic critical slope rule in 1+1 dimension |L7|. 
It has been shown p8[ that the same critical exponents 
occur in a wide class of critical systems, that includes 
the models of ricepile, depinning and earthquakes. The 
problem of an interface depinned from quenched defects 
in the anisotropic medium has been studied in |l9| ] . The 
exponents can be related to the directed percolation ex- 
ponents (see also p0|] ). It has been elucidated that due 
to anisotropy the threshold driving force depends on the 
interface slope, and thus may remain finite at the tran- 
sition when the velocity of the interface vanishes. The 
same effect is taken into account by the relevant nonlin- 
ear term \(dn h) 2 in the continuum equation of interface 
depinning pit]. In addition, the nonlinearity of the dy- 
namics in critical systems [001 leads to self-tuning to 
the dynamic critical state. 

Recently it has been shown jl6| that a self-tuning in the 
class of critical disordered systems can be described by 
self-consistent building of correlations in the transport 
equation with random quenched currents 



dh 



— = vndih + v ± dih - p{x)d\\h + n 
ot 11 



(2) 



Here p{x) is the random local velocity of transport with 
zero mean and correlations given by 
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{p(x)d ll hp(0)d ll h) d = 7 (^^j (d {l h) 2 , (3) 

and the anisotropy exponent £ is to be determined self- 
consistently as the system approaches the fixed point. 
The usual dynamic noise term 77 is assumed to have 
annealed correlations as (t?(%, t)rj(x't')) = 2DS^ d '(x — 
x')S(t — t'). It has been found in Ref. Em that a stable 
fixed point exists for finite strength of disorder 7 with 
the exponents that depend continuously on the range of 
correlations along parallel direction, parameter 5. In par- 
ticular, for long-range correlations at S = 1/3 the corre- 
sponding avalanche exponents are found r t = 1.89 and 

1/9 2 

t s = 1.56, with the emergent nonlinear term (dt\h) ■ 

x 11 



V. CONCLUSIONS 



eventually depinning of the wall occurs at critical driv- 
ing field. We suggest that in the zero driving rate the 
BN in this region belongs to the same universality class 
as the ricepile ||l7| and the anisotropic interface motion 
in 1 + 1 dimensions. For high disorder region, which 
is characterized by multidomain structure in the plane 
(see for instance fOg]), the criticality of BN is induced by 
the finite driving rates, which represents usually applied 
driving conditions in the experiments of BN. In princi- 
ple, a range of values of driving rates should exist for 
which a given sample exhibits self-similar Barkhausen 
noise. Both numerical simulations fL5| and experimen- 
tal observations P,[7|,p|JlO[ agree that scaling exponents 
decrease with increasing driving rate, with the linear de- 
pendence at large driving rates. Discussion of the finite 
driving rates in terms of the relevant fixed points of a 
transport equation has not yet received full theoretical 
treatment. 



Barkhausen avalanches in disordered ferromagnets ob- 
tained at sufficiently slow (which respects time scale sep- 
aration) increase of the driving field reflect fractal struc- 
ture of the dynamic critical states. The quantitative 
characterization of a given system in terms of sets of scal- 
ing exponents of the BN avalanche can be expected in 
several universality classes, which depend on: (i) type of 
domain walls; (ii) strength of pinning; and (Hi) applied 
driving conditions. The chemical composition, prepara- 
tion, sample form and annealing result in certain type of 
domain walls, which are thus characteristic for a given 
sample. Whereas, the driving rate can be controlled in 
the experimental setups. In ribbon geometry the effects 
of demagnetizing fields are minimized with longitudinal 
anisotropy [ ft2| . 

The scaling behavior studied in this work is expected to 
apply for the anisotropic samples with extended domain 
walls, which are pinned by nonmagnetic impurities. Such 
domain walls are found in ribbon samples with strong 
longitudinal anisotropy as for instance in fully annealed 
Fe78Bi3Sig ribbons in which domain walls tend to orient 
along longitudinal anisotropy axis 0,0. Similar struc- 
ture of domains was reported recently in Refs. ||[l(| for 
Fe2iCo64Bi 5 (see Fig. 1 in [@). Moreover, it has been 
suggested M that type of domains in these samples can 
be controlled by continuous variation of the tensile stress. 
Our results suggest (see also jll|) that self-similar 
Barkhausen noise in these systems occurs in a wide range 
of values of the disorder and driving field, which are 
marked as gray areas in the phase diagram in Fig. 2, 
excluding only the area with finite domain-wall velocity 
(upper left corner of the phase diagram). In the zero 
driving rate we find two distinct universal behaviors for 
low disorder and high disorder region, respectively, which 
compare well with the experimental results in the above 
samples. For low disorder self-similar BN is a conse- 
quence of the intermittent motion of a single domain wall 
(or well separated, noninteracting domain walls), before 



ACKNOWLEDGMENTS 

This work was supported by the Ministry of Science and 
Technology of the Republic of Slovenia. 



and A. Magni, J. Appl. Phys. 75, 
G. Bertotti, and A. Magni, Frac- 

Portteseil, J.Magn. Magn. Mat. 

and J. T. Markert, Phys. 



[1] L.V. Meisel and P.J. Cote, Phys. Rev. B 46, 10822 
(1992). 

[2] G. Bertotti, G. Durin, 
5490 (1994); G. Durin 
tals 3, 351 (1995). 

[3] O. Geoffroy and J.L. 
133, 1 (1994). 

[4] J. S. Urbach, R. C. Madison 
Rev. Lett. 75, 276 (1995). 

[5] G. Durin and A. Magni, in Fractal Reviews in the Natu- 
ral and Applied Sciences, edited by M. Novak (Chapman 
& Hall, London 1995), p. 35. 

[6] Dj. Spasojevic, S. Bukvic, S. Milosevic, and H. E. Stan- 
ley, Phys. Rev. E 54, 2531 (1996). 

[7] G. Durin, Proceedings of 14th International Conference 
on "Noise in Physical Systems and 1// Fluctuations", 
Eds. C. Claeys and E. Simoen, World Scientific, Singa- 
por, 1997. 

[8] G. Durin and S. Zapperi, e-print cond-mat/9808224 . 
[9] P. Bak, How Nature Works, (Springer- Verlag, New York 
1996). 

[10] S. Zapperi et al, Phys. Rev. B 58, 6353 (1998). 
[11] B. Tadic and U. Nowak, preprint. 

[12] JD. Livingston and WG. Morris, J. Appl. Phys. 57, 
3555 (1985). 

[13] E. Vives and A. Planes, Phys. Rev. E 50, 3839 (1994). 
[14] B. Tadic, Phys. Rev. Lett. 77, 3843 (1996). 
[15] B. Tadic, Philosophical Mag. B 77, 277 (1998). 
[16] B. Tadic, Phys. Rev. E 58, 168 (1998). 



4 



[17] A. Corral, Ph.D. Thesis, Barcelona University (1997). 
[18] M. Paczuski and S. Boettcher, Phys. Rev. Lett. 77, 111 
(1996). 

[19] L.-H. Tang, M. Kardar, and D. Dhar, Phys. Rev. Lett. 
74, 920 (1995). 

[20] B. Tadic and D. Dhar, Phys. Rev. Lett. 79, 1519 (1997). 
[21] M. Kardar, G. Parisi, and Y. C. Zhang, Phys. Rev. Lett. 
56, 889 (1986). 




■ I i . ■ ■ ■ - { ■ i ■ ■ in"" ■ i ■ i ■ ii. ^ vi ■ t: _!-T' 




FIG. 1. Cluster ol nipped (bright) and unflipped (gray) 
spins, with most recent cluster (dark gray) for L — 200 and 
/=0.3 in the presence of nonmagnetic defects (black points) 
of concentration c: Bottom: c=0.2; Top: c=0.05. 
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FIG. 2. Simulated Barkhausen noise signal vs time for one 
side of the hysteresis loop. Lattice size is L = 200 and dis- 
order values / = 0.3 and c = 0.2. Inset: Time series of the 
increments of the external field AH(t) vs t which are adjusted 
to the minimum local field in the system, corresponding to the 
signal in the main picture. Note that the main contribution 
to the distribution of avalanches comes from the central part 
of the hysteresis loop, which is reached after approximately 
500 steps. 
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FIG. 3. Phase diagram: external field H vs. concentration 
of nonmagnetic ions c for / = 0.3 and infinitely slow driving. 
Dashed line: Coercive field Hq(c) vs. c. Scale free Barkhausen 
avalanche (up to a finite cut-off) occur in the entire gray area, 
corresponding either to single-domain or multi-domain struc- 
ture. A single-domain depinning occurs along the boundary 
with finite velocity region. For a finite driving rate boundary 
between two scaling regions moves towards lower concentra- 
tions. 
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FIG. 4. Double logarithmic plot of the integrated distri- 
bution of avalanche durations P(t, c) for /=0.001 and var- 
ious concentrations c=0.1,0.125, 0.15, and 0.2, obtained by 
infinitely slow driving for L = 200. Inset: Critical exponents 
(A) Tt, and (o) r s vs. concentrations of nonmagnetic ions 
c. Full symbols: high disorder; Empty symbols: low disorder 
regime. Finite driving rate 5h = 0.01: (*) Tt, and (x) r s 
(from Ref. [14]). 



